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Abstract 

We consider the problem of covering hypersphere by a set of spherical hypercaps. This sort of prob- 
lem has mimcrous practical applications such as error correcting codes and reverse fc-ncarest neighbor 
problem. Using the reduction of non degenerated concave quadratic programming (QP) problem, we 
demonstrate that spherical coverage verification is NP hard. We propose a recursive algorithm based 
Z/1 on reducing the problem to several lower dimension subproblems. We test the performance of the 

I proposed algorithm on a number of generated constellations. We demonstrate that the proposed al- 

gorithm, in spite of its exponential worst-case complexity, is applicable in practice. In contrast, our 
■'"^ results indicate that spherical coverage verification using QP solvers that utilize heuristics, due to 

numerical instability, may produce false positives. 
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0\ 1 Introduction 
O 

We consider the problem of determining whether a hypersphere is completely covered by a set of hyper- 
spherical caps. An equivalent problem is whether the given set of hyperspherical cones centered in the 
origin, covers the whole space W^. A three dimensional version of the problem can be solved using the 
approach from [4], but a solution to a generalized problem in arbitrary dimensional space has not, to the 
5^ best of our knowledge, been proposed yet. The generalized problem of hyperspherical coverage by a set 
of hypercaps arises in areas such as coding theory [12] and multidimensional queries [11]. 

Covering problems are important in computational geometry and have been extensively studied re- 
cently. Elbassioni and Tiwary [8] considered the following problem: Given a set of polyhedral hyperconcs 
Ci,C2, . • • ,Cfc and a convex set D C W^, check whether cones cover the set D or not. They proved NP 
completeness in several cases and connected it to the problem of determining whether a union of the 
convex sets is convex. Also in [3] the authors considered covering a given set of points with a given 
polygon, whether in [5] authors considered covering the set of points with two disjoint disks and two 
disjoint squares. Papers [7] and [18] consider covering a sphere (ball) with other spheres (balls). 

Recently, development of algorithms for incremental density-based outlier detection [13, 16] have 
motivated the need to efficiently apply techniques for reverse fc-nearest neighbor search. The minimal 
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number of hypercaps that can completely cover a hypersphere is related to the theoretical upper bound for 
the number of reverse /c-nearest neighbors of a given point and hence the complexity of incremental outlier 
detection algorithms [15]. Also, the sets of hypercaps that can completely cover the hypersphere are basis 
for practical algorithms for reverse A;-nearest neighbors. In this paper, we consider hyper spherical coverage 
verification: for a given set of hypercaps on a hypersphere we determine whether the set completely covers 
the hypersphere. We are, however, not concerned how to determine the set of hypercaps that completely 
covers the hypersphere. 

Using the reduction of concave quadratic programming (QP) problem, we demonstrate that spherical 
coverage verification is NP hard. As a consequence, an algorithm with non-polynomial worst-case com- 
plexity may still be viable and practical. We provide a generalized recursive algorithm that can perform 
coverage verification task for arbitrary dimension d. The proposed algorithm is based on reducing the 
problem to several lower dimension subproblems. In addition, we provide a method that can identify a 
point on a hypersphere not covered by any hypercap, if such a point exists. 

We test the performance of the proposed algorithm on a number of generated constellations with 
diff'erent dimensionality. We demonstrate that the proposed algorithm, in spite of its exponential worst- 
case complexity, is applicable in practice, with acceptable average-case performance. In contrast, our 
results indicate that spherical coverage verification using heuristics-based QP solvers, may produce false 
positives and suffer from numerical instability. 

2 Spherical coverage verification 

In this section, we formally define a spherical coverage verification problem and demonstrate that the 
considered problem SphCovVer can be described as a system of non-linear equations and inequalities. 
Subsequently, we demonstrate that the problem at hand can be represented as quadratic programming 
problems with linear constraints. 

2.1 Problem formulation 

Suppose that we have n cones Ci, . . . , C„ in d-dimensional space W^, d <2 centered at point O = (0, . . . , 0) 
and defined by 

Ci = C{U-A) = |x G M-^ I > 9^ . 

Note that each cone Ci is defined by point ti G and real number — 1 < < 1. There holds 
cos ZxOti > 6i for each x e Ci and x O. 

For any two given points x,y G W^, with we denote usual scalar product as {x,y) = Yli=i^iyi 

and with ||a;|| we denote the Euclidian norm = a/ {x, x). 

Problem 1. (SphCovVer) Check if cones Ci cover the whole space W^. Equivalently, check if hypercaps 
Ki = CiD cover an unit hypersphere 5^(1) = {x € M*^ | = 1}. 

Without loss of generality, we can assume that all points ti belong to unit hypersphere £'^(1), i.e., 
which holds \\ti\\ = 1. Let ti = {tn, . . . , tid)- If x G Sd{l) then: 

d 

X e Ki <s=» ^ Xitij > 9i 
i=i 

Observe that point x on the unit hypersphere Sd-i{l) is not covered by any of the cones Ci, . . . , (7„ if 
and only if (x, ti) < 9i for alH = 1, . . . , n. Therefore cones Ci, . . . , C„ cover the unit hypersphere S^-ii^) 
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if and only if the following non-linear system of equations and inequities does not have a solution 

xl + xl + ... + xl = l 

tllXl + tl2X2 + . . . + tldXd < 01 

(1) 

tnlXl + tn2X2 + • • • + tndXd < Gn- 

2.2 Spherical coverage verification as quadratic optimization problem 

Denote by S the solution space of the linear system of inequalities obtained by dropping the first equation 
in (1). Obviously S is the convex set. Denote by S the closure of S i.e., the solution space of the 
inequalities from (1) when each < is replaced by <. The set S is also convex. Let f{x) = 
m = f{x) = min{/(x) | a; G 5} and M = f{x*) = max{/(x) | x G 5} (if 5 is unbounded, then 
M = +oo). In other words, x and x* (i.e., m and M) are solutions of the following QP problem: 

(min / max) f{x) = ||a;|p = x\+ x\ + . . . + x^ 

S.t. tiiXi + ti2X2 + . . . + tidXd < 6i 

(2) 

tnlXi + tn2X2 + ■ ■ • + tndXd < On- 

If S is unbounded (M = +oo) then let x* G 5 be any feasible point so that ||x|| > 1. We say that set 
S specified by constraints from (QP problem (2)) is degenerated if it is contained in some hyperplane T-L. 
Note that hyperplane % has to be of the form {x,ti) = 9i, i.e. there have to exist two constraints i and 
j from eq. (2) where t^k = —tjk, k = 1, . . . ,d and 9i = —6j. In such case, the system of equations (1) 
obviously has no solutions. 

The following lemma shows the connection between problem SphCovVer (i.e., the system (1)) and 
the QP problem (2): 

Lemma 1. The system of equations and inequities (1) has solutions if and only if M > 1, m < 1 and S 
is non-degenerated. 

Proof. 

(<^:) Let M > 1, m < 1 and S be non-degenerated. Assume that x* is a boundary point of S and 
consider a ball Bd{x*,p) where p < \\x*\\ — 1. Since 5 is a non-degenerated polytope, there exists an 
internal point G 5 fl Bd{x*,p). If x* is an internal point, we just set xl = x*. Note that p < \\x*\\ — 1 
implies ||x^|| > 1. 

The same way, we consider a boundary point x oi S and construct a new point xi G <S so that' 

||ii|| < 1 (here we take p < 1 — \\x\\). 

Since [ii,a;J] C 5 (5 is convex), function /(x) = is continuous on [xi,a:J] and f{xi) < 1 < f{x\), 
there exists point u G (xi, x\) C 5 so that f{u) = 1. In other words, u is the solution of system (1). 

(=>:) Let u be one solution of system (1). Hence u & S and f{u) = 1. Since S is an open set, there 
exists a ball Bd{u; p) C S. Denote by ui and U2 the intersection points of Bd{u] p) and line On, so 
that u G (Oui). Obviously holds f{ui) = (1 + p)"^ and /(«2) = (1 — pY- Now ui,n2 G Bd{u; p) C S 
and f{ui) > 1 > f{u2) directly implies M > f{ui) > 1 > f{u2) > m. Non-degeneracy of S follows 
immediately from the fact that (1) has solutions. □ 
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3 Spherical coverage verification is NP-hard 



In this section, we prove that the spherical coverage verification problem (SphCovVer) is an NP hard 
problem. First, we demonstrate that the concave non degenerated quadratic programming decision prob- 
lem (ConNDQPd) defined below can be polynomially reduced to SphCovVer. Then we demonstrate 
that the problem ConNDQPd is NP complete. 

Note that the variant of the problem ConNDQPd without the non-degeneracy assumption is con- 
sidered by Preund and Orlin in [10] (HB problem) where its NP completeness is proven. The degeneracy 
assumption makes the problem ConNDQPd considered here more restrictive than the one considered in 
[10], implying that we need a different proof of NP completeness. It will be given in the next subsection. 

3.1 Polynomial reduction of ConNDQPd to SphCovVer 

Consider the following concave quadratic programming (QP) problem. 
Problem 2. (ConNDQPd) Check whether exist xi, X2, . • . , a;,^ G M so that 

x-^ ~\~ X2 -|- • • . -|- x^ ^ c 
aiixi -I- ai2X2 . . . -I- aidXd < h 

(3) 



anlXi + an2X2 -|- • . . -|- ttndXd < K- 

where c, aij, 6i G M (c > 0) for i = 1, . . . ,n,j = 1, . . . ,d and the polytope specified by < constraints from 
(3) is non- degenerated. 

Also consider the following algorithm: 
Algorithm 1. QP-SphCovVer 

Input: An instance of the problem ConNDQPd, i.e., matrix A = [aij] G W^^'^ and vector b = [bi] € 

1. Normalize each constraint, i.e com,pute 

Uj = ."'^ ^ , 9i = —^ — i = l,...,n; j = l,...,d. (4) 

2. Solve the minimization problem (2) (a convex optimization problem) in polynomial time and denote 
its minimum by m. If m > 1, then output True. // the problem is infeasible, output False. 
Otherwise continue. 

3. Drop each constraint which satisfies 9i > 1. 

4. Form the instance of the problem SphCovVer from the remaining constraints and solve it. Output 
the complementary result. 

The following theorem proves the correctness of Algorithm QP- SphCovVer. 

Theorem 2. Algorithm QP-SphCovVer polynomially reduces the problem ConNDQPd to the problem 
SphCovVer. 

Proof. Note that by (4) wc form the equivalent problem of form (2) so that = 1 and c = 1. 

First assume that m > I. Then all the feasible points satisfy ||a;|p > 1 and hence the output of 
problem ConNDQPd is True. Assume that m = 1. By assumption, the feasible set of the problem 
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(3) is a non-degenerated polytope. Hence it cannot belong to the unit hypersphere (otherwise, it would 
reduce to the single point) and there is a feasible point x' so that > 1. It implies that the answer of 
ConNDQPd is True. 

Now assume that m < 1. Let x be the solution of minimization problem (2). Then the following 
holds 

\U\X\ +U2X2 + ... + tidXd\ = \iti,x)\ < \\x\\\\ti\\ = \\x\\ =m<l. (5) 

Equation (5) implies 6i > —1, since x is a feasible point. Without loss of generality, assume that 
01,62, . . . ,9p < 1 and ■ ■ ■ ,0n> ^- We can consider cones 

Ci = C{ti;9i), C2 = C{t2; O2), . . . ,Cp = C{tp; Op), 

since —l<Oi<l and \\ti\\ = 1 for all i = 1,2, . . . ,p. Also consider the corresponding system of first 
p+1 equations from (1): 

xl + xl + ... + xl = l 
hixi + ti2X2 + . . . + hdXd < Oi 

(6) 

tplXi + tp2X2 + . . . + tpdXd < Op. 

If cones Ci, C2, . . . ,Cp cover W^, it implies that the system (6) (and also (1)) does not have a solution. 
According to Lemma 1, there must hold M < 1, which implies that the answer of ConNDQPd is False. 

Assume that cones do not cover W'- and denote by x the solution of the system (6). Now, since the 
following relation holds for z = p + l,p + 2, . . . , n: 

tiixi + ti2X2 + . . . + tidXd = iU, x) < \\x\\ pill = 1 <0i, 

we conclude that x is also a solution of (1). According to Lemma 1, there holds M > 1 and the answer 
of ConNDQPd is True. This completes the proof. □ 



3.2 NP completeness of ConNDQPd and SphCovVer 

We prove that the problem ConNDQPd is NP hard by reducing it to the fc-clique decision problem. 
Recall that, for a given graph G = {V, E), set C 1/ is a clique, if for every u,v e CI holds {u, v} G E. 
In other words, a clique is every set CI of vertices, so that each two vertices from CI are adjacent. Clique 
CI is called k-clique, if it contains exactly k vertices. The A;-clique decision problem (see e.g., [6]) can be 
formulated as follows: 

Problem 3. (fc-Clique) Given a graph G = {V,E), check if there exists k-clique. 

It is known, ([6]) that the problem fc-Clique is NP complete. The following lemma demonstrates that it 
can be polynomially reduced to the problem ConNDQPd. 

Lemma 3. Problem /c-Clique can be polynomially reduced to problem ConNDQPd. 

Proof. For a given graph G = {V, E) with the vertex set V = {1,2,..., n}, consider the following 
instance of problem ConNDQPd: 

x\-\- x\ + . . . + x'^> n — e 

-l<Xi<l 

Xi^x^ < 0, V{i,j} ^E 
Xi + X2 + ■ . . + Xn > ^k — n 
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Here < e < 1 and its value will be determined later. If there exists a clique CI of length k in graph G, 
then by setting 

-1, i ^ CI 



Xi = 



we obtain one feasible solution {x\,X2, ■ ■ ■ , ) of the problem (7) satisfying (xl) + {X2) + . . . + (a;* ) = 
n > n — e. 

Now assume that there is no clique of length k in graph G. We show that the decision problem (7) 
does not have the solution. Consider the following auxiliary optimization problem 

max Xl + X2 + ■ ■ ■ + Xn 

s.t. — 1 < Xi < 1, 1 = 1,2, ...,n 
— 1—1 ) ) ) 

xf + X2 + ■ ■ ■ + x"^ > n — €. 

Let (xi,X2, ■ ■ ■ ,Xn) be an arbitrary feasible solution. It can be easily checked that > 1 — e (due to 
the quadratic condition in (8)) for i = 1,2,..., n implying that xi G [— 1, —\/l — e) U {s/l — e, 1]. Let 
Xi^,Xi^, . . . ,Xi^ G — e, 1] and Xq G [—1, — Vl — e) for q ^ {n,«2, • • • ,^«}- For arbitrary 1 < r < s < p 
there holds {i^, ig) G E, according to Xi^ + Xi^ > 2\Jl — e > 0. In other words, vertices 11,12, ... ,ip form 
a clique of length p in graph G. 

According to our assumption that there is no clique of length k in G, it must hold that p < k. 
Furthermore it holds that 

p 

f = Xl + X2 + . . . + Xn = '^Xij + ^ Xq <p - {n-p)y/l - e 

j=i qHii,i2,--4p} 

< p — {n — p){l — e) = 2p — n + e{n — p). 
Now by choosing e = 2/n and by using k > p and (9), wc obtain 

2 2(k — p) 
f < 2p — n-\ — (n — p) < 2p — n -\ (n — p) = 2k — n. 

n n — p 

According to the previous expression, each feasible solution of (8) satisfies f < 2k — n implying that the 
system (7) has no solutions. □ 

As a direct consequence of the Lemma 3 and the NP completeness of the problem fc-Clique, and 
since a verification of a solution for the eq. (3) is possible in polynomial time, the following corollary 
holds: 

Corollary 4. Problem ConNDQPd is NP complete. 

Now Corollary 4 and Theorem 2 directly imply: 
Theorem 5. Spherical coverage verification, i.e., the problem SphCovVer, is NP hard. 



4 Algorithms for spherical coverage verification 

The simplest method for spherical coverage verification is to apply a non-deterministic Monte-Carlo 
approach. The idea is to generate a large number N (for example N = 10^^) of pseudo-random points 
X distributed uniformly on the sphere. For each generated point x we check if there is hypcrcap Ki 
containing the point. If the answer for any point is negative, the algorithm outputs False. If the system 
of hypercaps K^, i = 1, . . . ,n covers the unit sphere, this method always outputs the correct answer True. 
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If there is no coverage, this method will output the correct answer False with probability that increases 
with A''. However, there is no guarantee that the algorithm will not return false positives (thus providing 
answer True for a non-covering system of hypercaps). Hence, in the subsequent subsections, we discuss 
algorithms for spherical coverage verification based on application of quadratic programming, and on the 
reduction of the given problem to lower dimensional subproblems. 

4.1 QP-based verification 

According to Section 2.1, cones Ci,i = 1, . . . ,n cover the space M*^ if an only if the system (1) has no 
solutions. According to Lemma 1, to check whether the system (1) has solutions, we need to solve the 
QP problems (2) and to determine whether the minimal value m and maximal value M of the objective 
function satisfy m < 1 < M. The following algorithm for solving problem SphCovVer arises from the 
aforementioned discussion: 

Algorithm 2. Cover-QP 

Input: Caps Ki, i = 1, . . . ,n defined by U G Sd{l) and 6i G (—1, 1). 

1. Return True if QP problems (2) are degenerated. Otherwise, solve both problems. 

2. If any of the problems is not feasible, or m < 1 < M does not hold, return True . Otherwise return 
False. 

To apply algorithm Cover-QP we need appropriate QP problem solvers. Note that the minimization 
problem can be solved in polynomial time, since it is convex. Since concave QP is NP complete (see 
Section 3.2), one of known heuristics can be applied ([17]). 

4.2 Recursive algorithm 

In this section, we describe our recursive algorithm for solving the problem SphCovVer. The main idea 
is to reduce the initial d-dimensional problem to several d — 1-dimensional problems. More precisely, 
using d-dimensional inversion, hypercaps are mapped into regions consisting of hyperspheres in d — 1 
dimensional plane and their exteriors/interiors depending on the position of the center of inversion w.r.t. 
a hypcrcap. As proven in the next subsection, a (i-dimcnsional sphere is covered by the hypercaps if and 
only if the resulting d—1 plane is completely covered by the d — 1 dimensional regions. In turn, this may 
be true if a boundary of each region (which is itself ad—1 dimensional hypersphere) is completely covered 
by d — 1 dimensional hypercaps defined by corresponding regions. Thus, the coverage of d-dimensional 
sphere reduces to coverage of d — 1 dimensional spheres. The following subsections provide the rationale 
for the algorithm and discuss the algorithm formally. 

4.2.1 Rationale of the algorithm 

We restate the well-known definition of the inversion in W^. 

Definition 1. Let c G be a given point, and let R be a positive real number. Inversion iPc,r{x) is a 
function il^c,R : — >■ M'^ that maps every point x to a point y so that: 

y = c+j. i^{x-c). (10) 

\\x — c|p 

Point c is called the center of inversion ipc,R o-iT'd R is the radius of inversion. 

Let us apply inversion tjj = ';/'(i,o,...,o),i caps Ki and unit hypersphere 5*^(1). It is well-known that 
an image of a hypersphere, by inversion whose center belongs to the hypersphere is a hyperplane. Thus, 
the image of the unit hypersphere is hyperplane xi = 1/2. 
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Denote by dX the boundary of a given set X, Particularly, we denote by Di = dKi the boundary of 
cap Ki. Also, for a given hypersphere S denote by intS* its interior and by extS* its exterior. Images of 
caps Ki, i = 1, . . . ,n consist of d — 1-dimensional hyperspheres, belonging to hyperplane xi = 1/2 and 
their exteriors or interiors, depending whether the center of inversion is outside or inside the cap. More 
precisely, the following Lemma holds: 

Lemma 6. Let c = (1, 0, . . . , 0) G M*^ and assume that c ^ Di for every i = 1, . . . ,n. Image of Di, by 
an inversion ipc,i is a d — 1-dimensional hypersphere Si with center Pi = {Pn, ■ ■ ■ ,Pid) O'^d radius ri. 
Moreover, the im,age of Ki is Ri = SiU extSj, if c E Ki and Ri = SiU intSi, if c ^ Ki. Values j3i and r^ 
are given by the following expressions: 

1 tij - 0^ 

= r = 2{9i- ti,y = 2{ei - ti,y ^^^^ 



Proof. Translate the coordinate system to the center of inversion, i.e., to point c = (1,0,..., 0). The 
equation describing cap Ki becomes {x' = x — c): 

tiix[ + ... + tidx'a>ei-tii. (12) 

Let tp{x) = y. By the involution property of inversion, we can conclude that ip{y) = x or in other words: 

1 , 

By replacing the last expression into (12) we obtain: 

d d 

Y.t,,y'^>{9.,~t,i)Y,yf- (13) 



yi = \- (14) 



Replacing y = y' + c finally yields: 

3 — 2 

Here condition 9i — > is equivalent to c G Ki. In (14) we denote: 

3 — ^ 

The last expression (15) can be further simplified using Yl'j=2'^lj — ^ ~ ^ii — h (H)- '-' 

In Fig. 1 we illustrate Lemma 6 for three spherical caps Di, i = 1,2,3 and their images. Observe 
that c & D3, hence D3 maps to the exterior of sphere S3. 

Due to Lemma 6, problem SphCovVer reduces to the following hyperplane cover verification (Hp- 
CovVer) problem: 

Problem 4. (HpCovVer) For given d — 1-dimensional hyperspheres Si, with center /5j and radius r^ 
belonging to the hyperplane xi = 1/2, and sets Ri so that Ri = SjUintSj or Ri = S'^ UextSj (i = 1, . . . ,n), 
check if sets Ri cover the whole hyperplane. 
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Figure 1: Three spherical caps and their images under the inversion ip for d = 3. 

To simphfy further discussion, if Ri = SiUintSi set Ri is called internal; otherwise we call it external. If 
all the sets Ri are internal, then we can immediately conclude that the answer to the problem HpCovVer 
is False. This holds obviously from the fact that all sets Ri are bounded and hence is also their union 



Suppose, in contrast, that at least one set Ri is external. If all sets R do not cover the whole 
hyperplane, there exists one hypersphere Si and point x £ Si so that it is uncovered by other sets Rj, 
j 7^ i. In other words, the following theorem holds: 

Theorem 7. Sets Ri, . . . ,Rn with different boundaries Si = dRi cover the whole hyperplane xi = 1/2 if 
and only if every hypersphere Si is covered by other sets Rj, j ^ i, i.e., 



Proof. 

(^:) If sets Ri,...,R n cover the whole hyperplane xi — 1/2, each hypersphere, Si as a subset of the 
hyperplane will be covered. 

(=^:) Denote an uncovered region of the hyperplane xi = 1/2 with Q. The boundary dQ consists of the 
union of spherical caps. Denote by A one of those caps and by Si the hypersphere which A belongs to. 
Since Si ^ Sj for i ^ j, interior intA cannot be covered by remaining hyperspheres Sj, j ^ i. In other 
words, there exists a point x £ mtA \ Uj^j ^j- Point x is covered by some set Rk and since x ^ Sk, there 
holds X £ intRk- Hence, there exists a ball 6) C intRk and holds B^_i{x; 5) n Q = 0. This is the 

contradiction with the fact that x G A is a, boundary point of Q. 



If Si = Sj for i ^ j then either Ri and Rj cover the hyperplane (if one of them is internal and the 
other is external) or Ri = Rj. In the second case, we can eliminate one of them and continue. 




(16) 



□ 
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According to Theorem 7 we need to check whether each hypersphere Si is covered by sets Rj, j ^ i. 
We distinguish the following cases, depending on whether the pairs of hyperspheres Si and Sj = dRj are 
disjoint: 

Case 1 . Hyperspheres Si and Sj have a nonempty intersection. In such a case Si H Rj is a hypercap Kj 
defined as: 

K] = \xeSi\i ^^^>m , ij^j, (17) 



^ 1 llx - /3i||||a;*-| 



where we define: 



2M— ' is internal ^ f d i(/3,- - ft), is internal , ^ no on 



2ridij 



, Rj is external ' ^ \-di/iPj - A)> is external 

(18) 

Observe that for the points Xj from eq. (18), x*- = {{Xj)2, ■ ■ ■ , ix^j)d) £ l^'' since (a;* )i = for every i 
and j. Case 1 is illustrated in the Fig. 2. 




Figure 2: Intersection of c? — 1-dimensional hyperspheres Si and Sj. 

Case 2. Hyperspheres Si and are disjoint. The equivalent condition is 5* ^ (—1,1), where is defined 
by eq. (18). In such case cither Si C Rj (Case 2a) or Si H Rj = 9 (Case 2b) holds. Which of these two 
sub cases holds can be determined e.g., by choosing the arbitrary point (for example Xj + (r^-, 0, . . . , 0)) 
on dSj and checking the inequality (14) for hypersphere Si. 

For fixed i, if for any j the condition Si C Rj (Case 2a) is satisfied, then eq. (16) holds and Si is 
covered. Therefore, the algorithm may continue with another value of i. Otherwise, it is sufficient to 
determine whether the sphere Si is covered by those hypercaps Kj defined by eq. (17) and corresponding 
to the pairs of spheres (Si,Sj) satisfying Case 1 (note that pairs {Si,Sj) satisfying Case 2b do not 
need be considered due to disjointness of Si and Rj). This is an instance of problem SphCovVer, for 
dimension d — 1. Hence, we reduce original problem SphCovVer to at most n — 1 equivalent problems 
of dimension d — 1. 
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4.2.2 Base case of the algorithm 

When d = 2, the inversion from Lemma 6 maps a 2D sphere (a circle) into a straight hne and 2D caps 
(arcs) degenerate into intervals. Hence, as the base case, we choose case d = 2 of problem HpCovVer. 
We omit the first coordinate (which is equal to 1/2) of each point from sets Si and Hence, we assume 
that there are given sets Ri so that Ri = [cj,di] {Ri is internal) or = M \ {ci,di) {Ri is external) for 
some real numbers Ci < di. Here q = /3j2 — Ti and di = j3i2 + rj. The problem is to check if sets Ri cover 
the whole real hne M, i.e., IJILi ^« ~ ^- Without loss of generality we can assume that . . . ,Rs be 
external and Rs+i, • • ■ , Rn internal. Let us define: 

c' = max Ci, d' — min d^. 

j=l,...,s i=l,...,s 

Obviously, if c' > c/', M is covered by the sets Ri. Otherwise, 

s 

{c!,d') = ^\[jRi, 
1=1 

and we need to check if the interval (c', d') is covered by [c^+i, ds+i], . . . , [c^, This can be performed 
by sorting segments [cs+i,ds+i], . . . ,[cn,dn] with respect to Ci and sequentially shortening the target 
interval (c', d'). Let us assume that the segments are sorted so that c^.+i < Cs+2 ■ ■ ■ ^ Cn- If c' < Cs+i 
then the interval (c', c^+i) is uncovered (since c^+i is minimal); hence M is also uncovered. Otherwise, we 
need to check if [c", d'], where c" = max{c', dg+i} is covered by segments [cs+2-, c?s+2], ■ ■ ■ ■, [cn, dn\. This 
leads to the following algorithm for solving the base case. 

Algorithm 3. Cover2-2D (d = 2 case of problem HpCovVer^ 

Input: Values j3i2, ri (vi > 0) and outj for i = 1, . . . ,n. We assume that out^ = True if Ri is external 
and otherwise outj = False. 

1. Let Ci = I3i2- ri and di = /3i2 + r, 

2. Reorder the sets Ri so that outj = True for i = 1, . . . , s and outj = False for ^ = s + 1, . . . , n. 

3. Sort sets Ri, i = 1 + s, . . . ,n so that Cs+i < . . . < c„. 
4- Let c' := maxj=i_..._s Cj and d' := minj=i^...^5 di 

5. For every i = s + 1, . . . ,n do the following: 

5.1. If c' > d' return True. Otherwise continue. 

5.2. If Ci > d then return False. Otherwise set c' := max{c!,di} and continue. 

6. If > d' return True, otherwise return False. 

Note that the complexity of the Algorithm Cover 2-2D is O(nlogn) if an asymptotically optimal 
sorting algorithm is used for intervals sorting. 

4.2.3 Algorithm outline 

Next, we formulate complete recursive Algorithm Cover for solving the general case of problem Sph- 
CovVer. 

Algorithm 4. Cover 

Input: Caps Ki, i = 1, . . . ,n defined by ti G Sd{l), {d > 2), and 9i G (—1, 1). 
1. If n = \ then return False. 
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2. Let c = (1, 0, . . . , 0) G M'^. Check if c G dKi for some i. In such a case, rotate the whole hypersphere 
in plane xiX2 by a small angle 6 so that c ^ dKi, i = 1, . . . ,n (condition of Lemma 6). 

3. Compute vectors Pi and values ri using eq. (11)- If tn > Oi (c G Ki) set out^ = True, otherwise 
set outj = False. // there holds out^ = False for every i = 1, . . . ,n then return False. 

4- If d = 2 apply Algorithm Cover2-2D for Pi2, ri and outj, i = 1, . . . ,n and return the obtained 
value. Otherwise continue. 

5. For every i = 1, . . . ,n do the following: 

5.1 Determine x^- and 6'j (relations (18) ), for every j ^ i, j = 1, . . . ,n. 

If 0* ^ (-1, 1), let A = x^j + (vi, 0, . . . , 0). If point A belongs to Rj (check the relation (14) ), 
set i = i + 1 and go to step 5. Otherwise continue. 

If 9] e (-1, 1), form cap K'-, eq. (17). If outj = True then set x) = -xj and 6) = -9]. 

5.2 Apply Algorithm Cover on the set of all formed caps Kj. 

Step 1 of Algorithm Cover implicitly covers the case when the number of caps n is smaller than the 
number of dimensions d. In step 2, wc introduce the rotation by a small angle S. We may set S to an 
arbitrary value, for example S = 0.01. If, after the rotation, point c is again on arc Cj, angle 6 needs 
to be changed. One possibility is to exponentially decrease it by setting 6 = pS, where < p < 1 (we 
used p = 0.9) and to repeat the same procedure until the point c is not on any arc Q. Step 3 performs 
inversion and checks whether the caps map into external or internal regions. If all regions are internal, 
as discussed earlier, the coverage is False. After this, we check whether the base case of recursion is 
achieved. Otherwise, Step 5 checks whether the conditions of the Theorem 16 are satisfied and, when 
needed, performs recursive calls of Cover. 

The worst-case time complexity of algorithm Cover is exponential in terms of d and polynomial in 
terms of n as shown by the following theorem: 

Theorem 8. The time complexity of Algorithm Cover is T{n,d) = O (n'''"^ logn) . 
Proof. 

Algorithm Cover, Step 5, reduces d-dimensional problem SphCovVer to at most n problems of 
size n — 1 and dimension d — 1. The complexity of reduce operations (including inversion, eq. (11), 

and checking conditions from Step 5), is 0{n^ • d). Hence, the following recursive relationship holds 
T{n, d) = 0{n'^ ■ d) + nT{n — l,d — 1). Since time complexity of the base case. Algorithm Cover2-2D, 
is ©(nlogn), this leads to the statement of the theorem. □ 

4.3 Localization of uncovered point 

When the solution to SphCovVer problem is False, it may be of interest to identify a point on a 
hypersphere not covered by any of the caps. We demonstrate how this could be accomplished using 
results of Algorithm Cover. 

The main idea of the proposed method is as follows. If unit hypersphere S'd(l), is not completely 
covered by caps Ki, . . . , then angles of the corresponding cones Ci can be slightly widened so that 
the resulting system of cones still does not cover the space. Further, as demonstrated in this section, it is 
possible to find a point on the boundary of the region covered by the enlarged cones, which corresponds 
to an uncovered point of the original problem. 

Namely, when the output of Algorithm Cover is False, we can determine a point u G 5*^(1) on the 
boundary of the covered region. Moreover, such a point belongs to an intersection of several boundaries 
Di = dKi, but u ^ intKi for every z = 1, . . . ,n. If such a point belongs to exactly I boundaries Di, we 
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call it (Z, d)-boundary point. We propose a method for computation of the {I, (i)-boundary point and then 
extend it to the computation of uncovered point. An example is illustrated in Fig. 3. The main idea is 
to identify a boundary point in a lower dimensional space (during recursive steps of Cover) and connect 
it with the boundary point of the original problem. 

Let the output of Algorithm Cover be False. Define recursively the sequence of points u"^ G >S'm(l) 
and v"^ e by 

where m = mo, mo + 1, . . . ,d, mo = maxjd — 1,2} and i is the index of uncovered hypersphere in step 5 
of m-th recursion call of Algorithm Cover. Also we denoted c"* = (1, 0, . . . , 0) and 



I 



d — 1, False is returned by Algorithm Cover 2-2D (step 4 of Cover), 
h — 1, False is returned at steps 1 or 3 of h-th recursion call of Cover. 



while /3™ and are corresponding values Pi and obtained from m-th recursion call of Algorithm 
Cover. In the first case {I = 1), initial point is equal to —1 or 1, depending of whether q or di is 
uncovered by sets Rj. Otherwise it is given by 



ci^h^i j -tf False is returned by Step 1 of Cover, 

c + , False is returned by Step 3 of Cover. 



Lemma 9, stated below, proves that u'^ is {I, d)-boundary point of the initial problem. 

Lemma 9. Every point w™, m = mo, mo + 1, . . . ,d, mo = min{d — l, 3} is a (m — d + l, m) -boundary point 
of the corresponding m,- dimensional SphCovVer problem in m-th recursion call of Algorithm Cover. 

Proof. We first identify a boundary point in the base case when the False answer of Algorithm 
Cover is detected. Then we prove by induction that each recursive call (step 5.2 of Cover) results in an 
additional boundary to which the point belongs. There are two possibilities for the base case of induction: 

1. Let the answer False be generated by Algorithm Cover2-2D. Define as follows 

1 I ~1) point Cj is uncovered 
u = < 

II, point di is uncovered 

It is not difficult to observe that = {1/2, ci) in the first and = (1/2, dj) in the second case. 
This point corresponds to point = ip^2 liv'^) G <S'2(1) (c^ = (IjO)) which is a boundary point of 
some arc Ki and uncovered by other arcs Kj, j / i. Hence, point ti^ is a required (1, 2)-boundary 
point for the case d = 2. 

2. Now let the answer False be generated by steps 1 or 3 of Algorithm Cover at recursion level h. 
Point u'^~f^+^ defined by eq. (19) does not belong to any hypercap (interior or boundary) and hence 
is uncovered. Therefore, u'^-''+^ is a {0,d — h-\- l)-boundary point. 

We prove the inductive step now. Consider the recursive call of Algorithm Cover on m-th level. If 
it is not otherwise stated, all notation corresponds to the SphCovVer problem being solved on m-th 
recursive level. fs 

Assume that the recursive call in step 5.2 returned False on the i-th. subproblem (i.e., covering of 
5m-i(l) by caps Kj, j = 1, . . . ,1). Also assume (by induction hypothesis) that returned point u"^~^ is 
a {m — 1 — d -\- l,m — l)-boundary point. Point v"^ = Pi -\- ri{0, u^"^) belongs to hypersphere Si. Since 
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jg — \ — d i^jji — l)-boundary point, does not belong to the interior of any cap Kf 

and hence does not belong to the interior of Rj (i.e., set Rj \ Sj) for any j. This directly implies that 
= ipc"^^i{v"^) is not contained in the interior of any hypercap Kj. 

Since point is, by assumption, contained in the boundary oim — 1 — d-\-l hypercaps Kj, there 
exist spheres 5^^, . . . , 5'i^_j_^^j so that G Si- for j = 1, . . . , Z — 1. Obviously, point n"* is contained 
within the boundary of / hypercaps Ki,Ki-^, . . . ■,Ki^_^_^^^. This completes the proof by induction that 
M*" is a (m — d + /, m)-boundary point. 

□ 




Figure 3: a) Spherical caps corresponding to cones Ci with 9 = \/3/2 (boundaries denoted by solid 
line) do not cover 53(1) and sperical caps (denoted by colored dashed lines) corresponding to cones Cf 
that also do not cover 53(1). b) Inversion '0(1,0,0),! of spherical cap boundaries; c) Inversion V'(i,o).i of 
2D sphere 5*2. Point v"^ is boundary point detected in Algorithm Cover2-2D corresponding to (1,2) 
boundary point and also to point G 5*2 and an uncovered point of 53(1). 

The following Lemma 10 formalizes the fact that given cones Ci, i = 1, . . . ,n which do not completely 
cover space M'^, can be enlarged (by increasing the central angle for a sufficiently small value) so that the 
resulting system of cones still does not cover the space. 

Lemma 10. // solution of SphCovVer problem is False for instance of cones Ci = C{ti;6i),i = 
l,...,n, then there exists sufficiently small value a so that' the solution 0/ SphCovVer problem for 
cones Cf = C{ti; 6i — a),i = 1, . . . ,n, is also False. 

Proof. Note that U = R'^\ Ur=i Ci 

is an open set and there exists at least one internal uncovered point 
X gU. Hence there exists a ball Bii{x; p) so that B^ix; p) C U. Let 7^ = arccos^j. All angles 7^ can be 
enlarged by value A7 = 2 arcsin(p/(2||x||)) and corresponding enlarged cones will not contain point x. 
Hence it is sufficient to set a = cos(A7). □ 

Note that (/, (i)-boundary point u'^'" of the enlarged coverage problem (from Lemma 10) is an internal 
uncovered point of the original problem. This holds since the distance between point u'^'°' and arbitrary 
cone Ci is at least 

,, / arccos(0j — a) — arccos(0j 

2 u ' sm ' 
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Hence, to find an uncovered point, it is sufficient to determine a (from Lemma 10), resolve the 
enlarged coverage problem and compute the boundary point (from Lemma 9). In practice, value a can 
be computed similarly as rotation angle 8 (see step 2 of Algorithm Cover). E.g., we can set a = 0.01 
and check if cones Cf,...,C^ cover the space W^. If the result is True, we can exponentially decrease a 
(e.g., by setting a = pa where < p < 1) and repeat the same procedure until the result is False. In 
our implementation, wc choose p = 0.9. 

Due to considerations above, the complete algorithm for computing an internal uncovered point can 
be formulated as follows: 

Algorithm 5. FindUncoveredPoint 

Input: Caps Ki, i = 1, . . . ,n defined by U G Sdil) and 9i G [—1, 1]. 

1. Apply Algorithm Cover with values ti and Oi- If the result is True, return True. Otherwise continue. 

2. Set a = 0.01. 

3. Set 9f = 9i — a. Apply Algorithm Cover with values ti and Of and compute {I, d) -boundary point 

4. If the result is True, set a = 0.9a and go to step 2. Otherwise return False and point u'^. 

5 Numerical examples 

In this section, we compare performance of proposed recursive algorithm Cover for spherical coverage 
verification, Section 4.2, with an algorithm based on quadratic programming Cover-QP, Section 4.1. 

We demonstrate that the application of Cover-QP could lead to false positives (coverage incorrectly 
verified) while the proposed recursive algorithm Cover docs not suffer from such a problem. Moreover, 
we demonstrate that the performance of Cover is satisfactory in practice, in spite of its worst-case 
exponential complexity. 

Algorithm Cover is implemented in programming language C. To implement Cover-QP, we utilized 
the programming packages Mathematica and Matlab. To test the influence of different non-convex QP 
solvers on the results of the algorithm, we also created an AMPL model [9] for the Algortihm Cover-QP 
and tested it using MINOS and FortMp solvers. 

Implementations are tested on several test examples. In the experiments, hypercaps are determined 
by constellations tj, i = 1 . . . , n of points and cones have a constant angle, i.e., Oi = 6. This stipulation 
comes from applications in methods for finding inverse fe-nearest neighbors. Namely, an algorithm for 
reverse A;-nearest neighbor problem from [1, 15] requires the covering constellation with minimal n and 
e = cos(^/6) = \/3/2. 

Test constellations are generated by the relaxation algorithm from [14]. This algorithm produces 
near-uniform placement of the points on d-dimensional hypersphere. It starts with a randomly generated 
set of initial points, where points interact through generalized electromagnetic interactions and each point 
has equal charge. The algorithm seeks the solution of the d-dimensional generalization of the Thomp- 
son's problem [2] , and searches iteratively for the equilibrium state (the state with minimal electrostatic 
energy) . 

5.1 Accuracy of algorithms 

Our experimental results indicate that Algorithm Cover-QP, which utilizes solvers for concave QP 
problems, can be numerically very unstable. As a consequence, the result is a potentially large number 
of false positives (an algorithm falsely indicates that a sphere is covered by caps). 

Consider constellation f our_D_85 obtained by relaxation algorithm for d = A and n = 85 (the whole 
constellation is given in the Appendix and can be found at tesla. cis .desu. edu/data/Constellations). 
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Mat lab implementation of the Algorithm Cover-QP returns True. We tested Mathematica implementa- 
tions for different working precisions (double precision, 20 digits, 50 digits and 100 digits). In all testings, 
result was True, but the corresponding optimal point and objective function values were different (see 
the following table). Results of testing AMPL model were similar. 





Objective 
function value 


Optimal point 


Xl 








Mat lab: double precision 


0.9756 


0.9477 


-0.0063 


-0.0004 


0.2314 


Mathematica: double precision 


0.935851 


0.172452 


0.898297 


-0.158036 


-0.272394 


Mathematica: 20 digits 


0.891469 


-0.390757 


0.136777 


0.77625 


0.342789 


Mathematica: 50 digits 


0.946319 


0.652345 


-0.476475 


0.236943 


0.487436 


Mathematica: 100 digits 


0.958959 


-0.901428 


0.034657 


-0.231813 


0.302405 


MINOS 


0.953615 


0.628285 


0.469643 


0.426044 


-0.39597 


FortMP 


0.947382 


0.037172 


0.869515 


0.0528789 


0.224396 



This example clearly shows the instability of Algorithm Cover-QP. One reason for such instability 
may be the fact that optimization heuristics embedded in the QP solver traps into the local maximum and 
do not achieve the global maximum. Such a conclusion is supported by the fact that the optimal solution 
found by the solver changes when simple variable transformation x = x' + c {c eW^ is a constant vector) 
is applied. For example, if we put c = (0.1,0.4,-0.7, —0.3) (this choice comes from the uncovered point 
shown in the next paragraph) all solvers return objective function value 1.000409053 and the optimal 
point (0.134335,0.4576476, -0.7915343, -0.3826162). Since the objective function value is greater than 
1, in such a case algorithm returns False. 

Note that Algorithm Cover-QP can have only false positives (for a non-covering constellation, pro- 
viding answer True), not false negatives. This can be explained by the fact that a point returned from 
the QP solver is always a feasible solution and the objective function value in this point cannot be larger 
than the global maximum. 

Unlike Cover-QP, Algorithm Cover is not based on heuristics and randomized initial conditions, but 
on relatively simple and numerically stable algebraic operations. During our experiments, we could not 
identify any case where Cover would return incorrect results. Particularly, on the same test example as 
above (the constellation f our_D_85) Cover correctly returns False, and detects the following uncovered 
point: 



X 



uncovered 



= (0.134309, 0.457496, -0.791181, -0.383002). 



Since maxi<j<85(a;"'^^°^^''^'^, ti) = 0.865901 < cos(7r/6) (see the Appendix for vectors U of the constellation 
f our_d_85) we can verify that 2;'^°^°™'"'''^ is indeed uncovered. 

Note that the result of the Monte-Carlo based approach method (see Section 4) on f our_D_85 with 
N = 10^° pseudo-random points was also incorrect (i.e.. True). This again demonstrates why Algorithm 
Cover seems to be the only reliable method for spherical coverage verification. 



5.2 Performance of recursive algorithm 

We test the working time of the Algorithm Cover. Implementation is compiled by a GNU C compiler and 
runs on the machine with AMD Phenom II CPU on 3.0 GHz and CentOS 5.2 (Linux) operating system. 
Testing is performed for d = 3, 4, 5 and 200 < n < 500. Results are shown in the Fig. 4. Execution 
times shown are obtained by averaging through 20 runs of the program on different constellations of the 
same dimensionality. All testing constellations corresponded to a sphere covered by caps (output of the 
Algorithm Cover was True). This restriction is added since the working time on constellations without 
coverage is considerably smaller due to the fact that Algorithm Cover does not complete recursion. 

As it can be seen from the graph, Algorithm Cover is practically applicable and does not reach 
its theoretically obtained complexity. This can be explained by the fact that not all d — 1-dimensional 
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Figure 4: Average running time of Algorithm Cover for d = 3,4, 5. 

hyperspheres Si intersect (Case 1 of subsection 4.1 holds). Practically, the number of hyperspheres 
intersecting Si is drastically smaller than n, which implies that the corresponding subproblems have a 
lower dimensionality. 

The applicability of Algorithm Cover made possible to determine the upper bound Mu((i) of minimal 
number M(d) of hypercaps with 9 = cos(7r/6), covering the unit hypersphere. In [15], it is proved that 
M(d) is at least Mi{d) = Q{2'^yfd), i.e., it is exponential as the function of d. In order to obtain Mu(d), 
we generate several constellations for fixed value n. If the covering constellation is found, n is decreased 
by 1 and the whole procedure is repeated. Otherwise we put 'M.u{d) := n. The following table shows the 
values of M.u{d) for d = 3,4, 5, 6: 



d 


3 


4 


5 


6 


M„(d) 


22 


81 


234 


715 



However, the closed-form expression for upper bound M.u{d) is still unknown and its existence is left 
as an open question. 

6 Conclusion 

Wc have considered the spherical coverage problem: given a set of hypercaps in c/-dimensional space, 
determine whether a d-dimensional hypersphere is completely covered by the hypercaps. We have demon- 
strated that the considered problem is NP hard by reducing concave quadratic programming (QP) prob- 
lem to it. We have discussed two algorithms to resolve the spherical coverage problem: the first method 
(Algorithm Cover-QP) is based on the utilization of quadratic programming. The second method (Algo- 
rithm Cover) is recursive and based on the reduction of the main problem on 0{n) problems of dimension 
d — 1. The recursive algorithm also provides a method to determine an uncovered point (if such a point 
exists). 

While the worst-case time complexity of the proposed recursive method is O (n''"-^ logn) , it is of 
practical interest, due to NP hardness of the considered problem (note also that the asymptotic complexity 
could be improved if, as the base case, we choose the case d = 3 of problem HpCovVer and utilize 
a method proposed in [4] to resolve it). However, numerical experiments indicate that the recursive 
algorithm almost never reaches maximal complexity and hence typically does not require prohibitive 
execution time. In contrast, Algorithm Cover-QP, with heuristics-based concave QP problem solvers, 
can be numerically unstable resulting in false positive detection (indicating false coverage of the sphere) 
and hence having limited practical application. 
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Our results indicate that the recursive algorithm may be the best algorithm for the problems having a 

relatively small dimension d. For the high values of d, where direct application of the recursive algorithm 
may be too time consuming, quadratic programming method still could be used as a presolve method, 
since it does not have false negatives. 

We conclude the paper with the following open problem: 

Problem 5. Develop an efficient algorithm for construction of the covering constellation of minimal 
size, for a given dimension d and angle 0. 
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A Constellation four_D_85 



tl 


t2 


<3 


ti 


0.911722 


0.083517 


-0.402106 


0.009974 


0.581301 


-0.242364 


0.278539 


-0.725096 


0.218664 


0.316284 


0.287611 


0.877172 


0.117213 


0.928431 


-0.32201 


0.143484 


-0.737931 


-0.158832 


-0.035717 


0.654946 


0.36765 


-0.515669 


-0.616088 


0.468352 


0.76061 


0.150628 


-0.006887 


0.631456 


-0.321449 


-0.24082 


0.870032 


0.285869 


0.344644 


0.293037 


0.546223 


-0.704976 


-0.412761 


0.540358 


0.013305 


-0.73312 


0.816541 


-0.164327 


0.467615 


0.295963 


-0.245525 


-0.585043 


-0.116853 


0.76406 


0.338649 


-0.178681 


-0.219111 


0.89743 


-0.69018 


-0.518921 


0.382654 


0.328555 


0.740476 


0.306078 


0.049102 


-0.596322 


-0.706023 


-0.187694 


0.659171 


-0.178316 


-0.055919 


-0.589895 


-0.805087 


-0.027052 


0.328685 


-0.068095 


0.878784 


0.339218 


0.943282 


-0.273977 


0.075663 


-0.171553 


-0.159834 


0.605736 


0.643967 


-0.43914 


0.059374 


0.684132 


0.716063 


0.125265 


0.562506 


0.316235 


0.756618 


-0.105413 


0.526948 


-0.692313 


-0.153046 


-0.468621 


0.585505 


-0.29157 


0.718579 


-0.236252 


-0.339912 


-0.865955 


-0.293451 


0.220156 
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tl 


t2 


^3 


t4 


0.504902 


0.689107 


-0.350765 


-0.383627 


-0.418404 


0.437219 


0.203124 


0.769752 


0.498145 


-0.17126 


-0.337765 


-0.780023 


-0.19152 


-0.885608 


-0.106055 


-0.409599 


-0.216553 


0.126561 


-0.29997 


0.920383 


0.060341 


0.820685 


0.197331 


0.532819 


-0.736942 


-0.41477 


-0.460301 


0.270194 


-0.547833 


-0.450474 


-0.569484 


-0.415501 


0.204308 


-0.777431 


0.433251 


-0.407619 


-0.983534 


-0.127577 


0.121635 


0.039876 


-0.106398 


0.780104 


-0.465982 


-0.403706 


-0.824089 


0.406231 


0.311938 


-0.241967 


0.925197 


0.302318 


0.226185 


-0.038143 


0.582566 


0.128159 


-0.589217 


0.544991 


0.622845 


-0.704481 


0.326492 


0.095782 


-0.836308 


0.452032 


-0.011796 


0.310027 


-0.309656 


0.548631 


-0.775807 


0.035212 


-0.225897 


-0.143505 


0.308569 


0.912777 


0.656179 


-0.438453 


-0.607877 


-0.087605 


-0.12433 


-0.638736 


0.476572 


0.591132 


-0. 425250 


().29«7ir) 


0.85 1262 


-0.012046 


0.480447 


0.819603 


0.259428 


-0.173544 


-0.514513 


0.096686 


0.503361 


-0.687428 


0.267199 


0.5784 


-0.305584 


0.707585 


-0.018035 


-0.020976 


-0.880919 


-0.472457 


-0.736443 


-0.653951 


-0.063084 


-0.161308 


0.033759 


-0.030823 


0.959959 


-0.276386 


-0.471801 


0.343863 


-0.601331 


-0.545493 


-0.349645 


0.727398 


-0.321563 


0.495215 


0.004942 


-0.96662 


0.232848 


0.106783 


0.111029 


-0.625644 


0.768139 


0.078766 


0.060656 


-0.490269 


-0.013342 


-0.869356 


0.588121 


0.537068 


0.416267 


0.438626 


0.338081 


-0.727789 


0.016177 


0.596458 


-0.342943 


-0.152519 


-0.336818 


-0.863529 


-0.868041 


0.131346 


-0.46866 


-0.098029 


0.084497 


-0.533808 


-0.58593 


-0.603817 


0.555933 


0.174623 


-0.695937 


-0.419663 


0.409128 


-0.265222 


0.414691 


0.768312 


-0.3i)cS()58 


-().719()19 


0.52()1()8 


-0.2i()135 


-0.230651 


-0.326573 


-0.707467 


0.582787 


-0.629986 


0.173969 


-0.5524 


0.517404 


-0.389677 


-0.043493 


-0.919642 


0.0228 


-0.013034 


0.298631 


-0.792652 


0.531369 


0.186482 


0.689421 


0.091131 


-0.693987 


0.814852 


-0.386301 


-0.176673 


0.39443 


0.275837 


-0.088184 


-0.946315 


0.143616 


-0.497203 


-0.491744 


0.211614 


-0.682786 


-0.11172 


0.959906 


0.17654 


-0.186902 


-0.028234 


0.117661 


0.115501 


-0.98591 


0.083739 


0.314023 


-0.444543 


-0.834721 


-0.805373 


-0.010126 


-0.060896 


-0.589545 


-0.125314 


0.270743 


0.737091 


0.606376 


0.690523 


0.637669 


-0.183595 


0.287835 


-0.481723 


0.736347 


0.436683 


0.1872 


-0.719834 


0.10622 


0.559014 


0.397567 


0.38002 


0.543697 


-0.745668 


0.062905 


-0.653788 


0.707026 


-0.192191 


-0.189048 
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tl 


t2 


^3 


t4 


0.013797 


-0.290263 


0.625477 


-0.72411 


0.303814 


-0.890275 


-0.328701 


0.084045 



Pseudocode of Algorithm Cover 

we give the complete pseudo-code description of Algorithm Cover. 
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Algorithm 1 Cover(;;; J: X; (9) - Splic-rical (-(n'orMgc \'('rilicali 



9 
10 
11 
12 
13 
14: 
15 
16 
17; 
18 

19: 

20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34: 
35 
36 
37: 
38 
39 
40 
41 
42 
43 
44: 
45 
46 
47: 
48 
49 



if d = 2 then 

return Cover2(n; X; 6') 
end if 

S := (1,0,..., 0) eW^ 

if {S, ti) = 9i for some i = 1,2, . . . ,n then 

X := Rot(n; d; X; 6) (rotates by small angle 6). 
end if 

ri := 2(g._f .'^) for i = 1, ■ ■ ■ ,n (radii of the spheres in d — 1-dimensional hyperplane) 

j3ij := 2(e'^-ti) i = 1, . . . ,n and j = 2,. . . ,d (centers of the spheres in d— 1-dimensional hyperplane) 

outj := {{S, ti) < di) (True if S belongs to the interior of cap Ki and False otherwise) 

for i := 1 to n do 

A; := 

Coverj := False 
for j := 1 to n do 
if i / J then 

k:=k + l 

tli_i ■.= Pji-l3ii for Z = 2,3,...,d 

d'- — Wf^W 

•— ll^fcll ^ ^ 

dij 

a -l<el<l then 
if outj then 

th.= -tlifovl = l,2,...,d-l 
end if 
else 

if (s < rj and not outj) or (s > rj and outj) then 

Cover, := True (j-th hypersphere completely covers i-th) 

break 
else 

k := k — 1 {j-th. hypersphere is disjoint from i-th) 
end if 
end if 
end if 
end for 

if not Coverj then 
ii k = then 

return False 
else 

Xi := [(ti)r...(ti)r]r g^fcxd-i 

0' ■■= i0l,---,el) 

Covers := Cover(A;; d - 1; X^; 9^) 
if not Coverj then 

return False 
end if 
end if 
end if 

end for 22 
return True 
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